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Abstract. Using existing models for the evolution of clustering we show that the 
first sources are likely to be highly biased and thus strongly clustered. We compute 
intensity fluctuations for clustered sources in a uniform IGM. The strong fluctuations 
imply the universe was reheated and reionised in a highly inhomogeneous fashion. 



INTRODUCTION 

Fluctuations in the ionizing intensity during reionization result in temperature 
and ionization fraction fluctuations. This has important implications for the reheat- 
ing and reionization history of the IGM. There are two contributing components: 

• IGM density fluctuations (variable optical depth) 

• Spatial distribution of ionizing sources 

The second point is usually addressed by computing Poisson fluctuations in the 
ionizing intensity ( [1], [2]). Here we consider a hitherto neglected effect, the clus- 
tering of ionizing sources, and show that it is the dominant source of intensity 
fluctuations. 



MODELLING THE IONIZING SOURCES 

We work within the framework of the concordance 
cosmology of [3] :(n m ,n A ,n b ,h,a 8h -i,n) = (0.35,0.65,0.04,0.65,0.87,0.96). We 
parametrise our model of reionization with 3 relevant lengthscales: 

•Mean separation n^ 1 ^ 3 The comoving number density of collapsed objects above 
a critical mass M* is given by Press-Schechter theory: n{z) = At hps(M, z)dM, 
where the creation rate of dark matter halos is given in [4]. The critical mass 
M* is given by the Jeans mass after reheating, Mj 

eans r * J 

10 9 (1 + z/10)- 3 / 2 M Q . 

Before reheating, it is given by the virial mass necessary for efficient atomic cooling 
M coo i ~ 10 8 (1 + z/10y 3/2 M Q . The source lifetime is assumed to be At ~ 10 7 yrs, 
the characteristic timescale for both starbursts and quasars. 

• Attenuation length r atten This is modelled as mm(r strom , r Hfe , r Lyman ), where 
r strom is the Stromgren radius, r^/e is the radius of the volume ionized by a source 




FIGURE 1. The evolution of various lengthscales with redshift. The increasing bias factor with 
rcdshift counteracts the fall of the linear growth factor, causing r Q to remain almost constant. 
The downturn of r atte „ at low redshift is caused by the increased clumping of the gas, which 
shortens the recombination time. 

during its lifetime, and r Lyman is the mean free path of photons bounded by the 
covering factor of Lyman limit systems. We assume values for the clumping of the 
IGM computed in [5]. The covering factor of Lyman limit systems is computed by 
assuming a mini-halo model (e.g., [6]), and computing their abundance via Press- 
Schechter theory. 

• Correlation length r Q The evolution of the matter correlation function with 
redshift is given by integrating the power spectrum. We assume a linear bias model 
6ih( r ) — fr 2 £mm( r )- Given the linear bias factor b(M) as in [7], we then compute 
6(M, z) by assuming linear theory and thus compute a number weighted bias factor 
at each epoch b(z): 

[Z) J£n(M,z)dM [i) 

The high bias means that the correlation length of halos is large, even though 
the correlation length of matter falls with redshift. For example, b ~ 7 at z=12, 
which leads to r)j al0 ~ S.SMpc/i" 1 in comoving coordinates, only somewhat weaker 
than the clustering of objects seen today. We also calculate the evolution of the 3 
and 4 point correlation functions, assuming standard hierarchical scaling relations, 
£iv( r ) = Sat^ -1 ^). Note that the Sn, although scale invariant, are a function of 
bias as well. 

To estimate the fluctuations in the radiation field, we conduct a Monte-Carlo sim- 
ulation, in the spirit of [1]. We generate a mock catalog of ionizing sources in 
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FIGURE 2. The probability distribution function of counts in cells at z ~ 12 for clustered and 
Poisson distributed sources, for a sphere of radius r atte „ ~ 4 Mpc (comoving). The clustered case 
has a high N tail, where many low luminosity sources can mimic an ultra luminous source. Note 
also the greatly increased probability for N=0 (voids) in the clustered case. 

a box 150 Mpc comoving size at z reion ~ 12. A fractal model [8] is used to lay 
down sources with specified 2, 3 & 4 point correlation functions, with appropriate 
modelling of the bias factors. This has the chief advantage of speed and conve- 
nience over using N-body simulations. It also does not suffer from mass resolutions 
effects in identifying halos. We verify that the 2 point correlation function and 
the 2nd, 3rd and 4th moments of counts in cells (which are directly related to the 
volume-averaged correlation functions, with corrections for shot noise) are correctly 
recovered. Luminosity is then randomly assigned to each source, assuming a Press- 
Schechter mass function, and a constant emissivity per unit mass, as given by [9] 
for their mini-quasar model. The ionizing flux at random points in the box is then 
computed, assuming each source ionizes out to a radius r atten . We thus compute 
the distribution function of the ionizing flux, as well as the two point correlation 



function OjW =< (J(0)~ < J >)(A r )- < J >) > / < J > 2 - 

CONCLUSIONS & OBSERVATIONAL IMPLICATIONS 



Clustering of ionizing sources produces a large dispersion in the luminosity of 
random cells. Clustered sources collectively mimic an ultra-luminous source. In 
addition, large voids exist where no ionizing sources are present. These are ionized 
by the percolation of ionizing radiation. Thus, large intensity fluctuations arise in 
the radiation field, which show coherence on large scales. 
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FIGURE 3. The two-point correlation function of ionizing intensity at z ~ 12 as a function of 
separation, assuming a sharp cutoff in ionizing flux at r atte „. In the clustered case, the radiation 
field displays correlations over large lcngthscales. 

This has important observational implications. Emission line searches (e.g. in 
Lya and Ha) at high redshift have a stronger signal due to the higher intensity 
fluctuations. Furthermore, CMB fluctuations due to Thomson scattering off elec- 
trons in moving ionized patches [10] are enhanced, both due to the larger patch 
sizes and the spatial correlations of the ionized patches. Finally, uneven reheating 
provides a possible biasing mechanism. 
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